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Glossary 

AKNS method A method introduced by Ablowitz, Kaup, Newell, and Segur in 1973 that identifies the 
nonlinear partial differential equation (NPDE) associated with a given first-order system of linear ordinary 
differential equations (LODEs) so that the initial value problem (IVP) for that NPDE can be solved by the 
inverse scattering transform (1ST) method. 

Direct scattering problem The problem of determining the scattering data corresponding to a given 
potential in a differential equation. 

Integrability A NPDE is said to be integrable if its IVP can be solved via an 1ST. 

Inverse scattering problem The problem of determining the potential that corresponds to a given set of 
scattering data in a differential equation. 

Inverse scattering transform A method introduced in 1967 by Gardner, Greene, Kruskal, and Miura that 
yields a solution to the IVP for a NPDE with the help of the solutions to the direct and inverse scattering 
problems for an associated LODE. 

Lelx method A method introduced by Lax in 1968 that determines the integrable NPDE associated with a 
given LODE so that the IVP for that NPDE can be solved with the help of an 1ST. 

Scattering data The scattering data associated with a LODE usually consists of a reflection coefficient 
which is a function of the spectral parameter A, a finite number of constants Xj that correspond to the 
poles of the transmission coefficient in the upper half complex plane, and the bound-state norming constants 
whose number for each bound-state pole Xj is the same as the order of that pole. It is desirable that the 
potential in the LODE is uniquely determined by the corresponding scattering data and vice versa. 

Soliton The part of a solution to an integrable NPDE due to a pole of the transmission coefficient in the 
upper half complex plane. The term soliton was introduced by Zabusky and Kruskal in 1965 to denote a 
solitary wave pulse with a particle-like behavior in the solution to the Korteweg-de Vries (KdV) equation. 
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Time evolution of the scattering data The evolvement of the scattering data from its initial value 
S{X, 0) at f = to its value ^(A, t) at a later time t. 

I. Definition of the Subject and Its Importance 

A general theory to solve NPDEs does not seem to exist. However, there arc certain NPDEs, usually 
first order in time, for which the corresponding IVPs can be solved by the 1ST method. Such NPDEs are 
sometimes referred to as integrable evolution equations. Some exact solutions to such equations may be 
available in terms of elementary functions, and such solutions are important to understand nonlinearity 
better and they may also be useful in testing accuracy of numerical methods to solve such NPDEs. 

Certain special solutions to some of such NPDEs exhibit particle-like behaviors. A single-soliton solution 

is usually a localized disturbance that retains its shape but only changes its location in time. A multi-soliton 
solution consists of several solitons that interact nonlinearly when they are close to each other but come out 
of such interactions unchanged in shape except for a phase shift. 

Integrable NPDEs have important physical applications. For example, the KdV equation is used to 
describe [14,23] surface water waves in long, narrow, shallow canals; it also arises [23] in the description 
of hydromagnetic waves in a cold plasma, and ion-acoustic waves in anharmonic crystals. The nonlinear 
Schrodinger (NLS) equation arises in modeling [24] electromagnetic waves in optical fibers as well as surface 
waves in deep waters. The sine-Gordon equation is helpful [1] in analyzing the magnetic field in a Josephson 
junction (gap between two superconductors). 

II. Introduction 

The first observation of a soliton was made in 1834 by the Scottish engineer John Scott Russell at the 
Union Canal between Edinburgh and Glasgow. Russell reported [21] his observation to the British Associa- 
tion of the Advancement of Science in September 1844, but he did not seem to be successful in convincing 
the scientific community. For example, his contemporary George Airy, the influential mathematician of the 
time, did not believe in the existence of solitary water waves [1]. 

The Dutch mathematician Korteweg and his doctoral student de Vries published [14] a paper in 1895 
based on de Vries' Ph.D. dissertation, in which surface waves in shallow, narrow canals were modeled by 
what is now known as the KdV equation. The importance of this paper was not understood until 1965 even 
though it contained as a special solution what is now known as the one-soliton solution. 

Enrico Fermi in his summer visits to the Los Alamos National Laboratory, together with J. Pasta and S. 
Ulam, used the computer named Maniac I to computationally analyze a one-dimensional dynamical system of 
64 particles in which adjacent particles were joined by springs where the forces also included some nonlinear 
terms. Their main goal was to determine the rate of approach to the equipartition of energy among different 
modes of the system. Contrary to their expectations there was little tendency towards the equipartition of 
energy but instead the almost ongoing recurrence to the initial state, which was puzzling. After Fermi died 
in November 1954, Pasta and Ulam completed their last few computational examples and finished writing 
a preprint [11], which was never published as a journal article. This preprint appears in Fermi's Collected 
Papers [10] and is also available on the internet [25]. 

In 1965 Zabusky and Kruskal explained [23] the Fermi-Pasta-Ulam puzzle in terms of solitary wave 
solutions to the KdV equation. In their numerical analysis they observed "solitary-wave pulses," named 
such pulses "solitons" because of their particle-like behavior, and noted that such pulses interact with each 
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other nonlinearly but come out of interactions unaffected in size or shape except for some phase shifts. Such 
unusual interactions among soHtons generated a lot of excitement, but at that time no one knew how to solve 
the IVP for the KdV equation, except numerically. In 1967 Gardner, Greene, Kruskal, and Miura presented 
[12] a method, now known as the 1ST, to solve that IVP, assuming that the initial profile u{x, 0) decays to 
zero sufficiently rapidly as x ^ ±oo. They showed that the integrable NPDE, i.e. the KdV equation, 

Ut - 6uUx + Uxxx = 0, (2.1) 
is associated with a LODE, i.e. the 1-D Schrodinger equation 

-0+«(a;,i)V = fcV, (2.2) 

and that the solution u{x,t) to (2.1) can be recovered from the initial profile u(.t,0) as explained in the 
diagram given in Section III. They also explained that soliton solutions to the KdV equation correspond to 
a zero reflection coefflcient in the associated scattering data. Note that the subscripts x and t in (2.1) and 
throughout denote the partial derivatives with respect to those variables. 

In 1972 Zakharov and Shabat showed [24] that the 1ST method is applicable also to the IVP for the 
NLS equation 

iut + Uxx + ^1^"^^- = 0, (2.3) 



where i denotes the imaginary number y— 1. They proved that the associated LODE is the first-order linear 
system 

— = —tXt, + U(X, t) 7], 

dx 

(2.4) 

— =iXr]-u{x,t)^, 

where A is the spectral parameter and an overline denotes complex conjugation. The system (2.4) is now 
known as the Zakharov-Shabat system. 

Soon afterwards, again in 1972 Wadati showed in a one-page publication [22] that the IVP for the 
modified Korteweg-de Vries (mKdV) equation 

ut + Gu'^Ux + Uxxx = 0, (2.5) 

can be solved with the help of the inverse scattering problem for the linear system 

d£ , _ , . 

— = -tX^ + u[x,t)r], 

(2.6) 

drj / X ^ 

— = iXri — ulx, t) t. 

- dx 

Next, in 1973 Ablowitz, Kaup, Newell, and Segur showed [2,3] that the IVP for the sine-Gordon equation 

Uxt = sinu, 

can be solved in the same way by exploiting the inverse scattering problem associated with the linear system 

— = -iX^- -Ux{x,t)i], 

dr] 1 . . ^ 

■^=iXv+^Ux{x,t)^. 
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Since then, many other NPDEs have been discovered to be solvable by the 1ST method. 

Our review is organized as follows. In the next section we explain the idea behind the 1ST. Given a 
LODE known to be associated with an integrablc NPDE, there are two primary methods enabling us to 
determine the corresponding NPDE. We review those two methods, the Lax method and the AKNS method, 
in Section IV and in Section V, respectively. In Section VI we introduce the scattering data associated with 
a LODE containing a spectral parameter and a potential, and we illustrate it for the Schrodinger equation 
and for the Zakharov-Shabat system. In Section VII we explain the time evolution of the scattering data 
and indicate how the scattering data sets evolve for those two LODEs. In Section VIII wc summarize the 
Marchenko method to solve the inverse scattering problem for the Schrodinger equation and that for the 
Zakharov-Shabat system, and we outline how the solutions to the IVPs for the KdV equation and the NLS 
equation are obtained with the help of the 1ST. In Section IX we present soliton solutions to the KdV and 
NLS equations. A brief conclusion is provided in Section X. 

III. Inverse Scattering Transform 

Certain NPDEs arc classified as integrablc in the sense that their corresponding IVPs can be solved with the 
help of an 1ST. The idea behind the 1ST method is as follows: Each integrablc NPDE is associated with a 
LODE (or a system of LODEs) containing a parameter A (usually known as the spectral parameter), and the 
solution u{x,t) to the NPDE appears as a coefficient (usually known as the potential) in the corresponding 
LODE. In the NPDE the quantities x and t appear as independent variables (usually known as the spatial 
and temporal coordinates, respectively), and in the LODE x is an independent variable and A and t appear 
as parameters. It is usually the case that u(x,t) vanishes at each fixed t as x becomes infinite so that a 
scattering scenario can be created for the related LODE, in which the potential u{x,t) can uniquely be 
associated with some scattering data ^(A, t). The problem of determining ^(A, t) for all A values from u{x, t) 
given for all x values is known as the direct scattering problem for the LODE. On the other hand, the 
problem of determining u{x,t) from S{X^t) is known as the inverse scattering problem for that LODE. 

The 1ST method for an integrablc NPDE can be explained with the help of the diagram 

f direct scattering for LODE at t=0 

u{x,Q) > 

solution to NPDE 

u{x,t) < 

inverse scattering for LODE at time t 

In order to solve the IVP for the NPDE, i.e. in order to determine u{x, t) from u{x, 0), one needs to perform 
the following three steps: 

(i) Solve the corresponding direct scattering problem for the associated LODE at t = 0, i.e. determine the 
initial scattering data 5(A, 0) from the initial potential u{x, 0). 

(ii) Time evolve the scattering data from its initial value S{\Qi) to its value S{X,t) at time t. Such an 
evolution is usually a simple one and is particular to each integrablc NPDE. 

(iii) Solve the corresponding inverse scattering problem for the associated LODE at fixed t, i.e. determine 
the potential u{x,t) from the scattering data S{X,t). 

It is amazing that the resulting u{x,t) satisfies the integrablc NPDE and that the limiting value of u{x,t) 
as t ^ agrees with the initial profile u{x, 0). 



5(A,0) 

time evolution of scattering data 

S{X,t) 
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IV. The Lax Method 

In 1968 Peter Lax introduced [15] a method yielding an integrable NPDE corresponding to a given 
LODE. The basic idea behind the Lax method is the fohowing. Given a hnear differential operator C 
appearing in the spectral problem Lil) = A^', find an operator A (the operators A and £ are said to form a 

Lax pair) such that: 

(i) The spectral parameter A does not change in time, i.e. At = 0. 

(ii) The quantity -i/'t — A'^l) remains a solution to the same linear problem CaI) = Xip. 

(iii) The quantity £t + CA — AC is a multiplication operator, i.e. it is not a difi^erential operator. 
From condition (ii) we get 

C{^t-AiP) = \{'^t-AiP), (4.1) 
and with the help of = X%1} and A* = 0, from (4.1) we obtain 

£Vt - = X'il)t-A (AV-) = dt - ACip = dt {Cip) - ACip = Cttp + C^Jt - AC^, (4.2) 

where dt denotes the partial differential operator with respect to t. After canceling the term Cipt on the left 
and right hand sides of (4.2), we get 

{Ct + CA-AC)i) = 0, 

which, because of (iii), yields 

Ci + CA- AC = 0. (4.3) 

Note that (4.3) is an evolution equation containing a first-order time derivative, and it is the desired integrable 
NPDE. The equation (4.3) is often called a compatibility condition. 

Having outlined the Lax method, let us now illustrate it to derive the KdV equation (2.1) from the 
Schrodinger equation (2.2). For this purpose, we write the Schrodinger equation as Cij) = Xip with A := fc^ 
and 

£:=-dl+u{x,t), (4.4) 

where the notation := is used to indicate a definition so that the quantity on the left should be understood 
as the quantity on the right hand side. Given the linear differential operator £ defined as in (4.4), let us try 
to determine the associated operator A by assuming that it has the form 

A = a3dl + a2dl + aid^ + ao, (4.5) 

where the coefficients aj with j = 0, 1, 2, 3 may depend on x and t, but not on the spectral parameter A. 
Note that Ct = Ut- Using (4.4) and (4.5) in (4.3), we obtain 

{)dl + {)dt + {)dl + {)dl + {)d, + {)=0, (4.6) 

where, because of (iii), each coefficient denoted by ( ) must vanish. The coefficient of vanishes automati- 
cally. Setting the coefficients of 9| to zero for j = 4, 3, 2, 1, we obtain 

3 3 

as = ci, (22 = C2, ai = C3 - -ciu, ao = C4 - -ciu^ - C2U, 

with ci, C2, C3, and C4 denoting arbitrary constants. Choosing Ci = —4 and C3 = in the last coefficient in 
(4.6) and setting that coefficient to zero, we get the KdV equation (2.1). Moreover, by letting C2 = 04 = 0, 
we obtain the operator A as 

A = -4dl + 6ud^ + 3u^. (4.7) 
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For the Zakharov-Shabat system (2.4), we proceed in a similar way. Let us write it as £V = -^V'j where 
the Hnear differential operator £ is defined via 





"1 " 




u{x,t) 


£ :=i 




dx - i 




-1 




u{x,t) 



Then, the operator A is obtained as 





"1 


■ 




"0 


u 




Ux 


A = 2i 





-1 


dl-2i 


u 





dx - i 


Ux 



(4.8) 



and the compatibility condition (4.3) gives us the NLS equation (2.3). 

For the first-order system (2.6), by writing it as £tp = Xtp, where the linear operator £ is defined by 





"1 " 




u{x,t) 


C :=i 




dx - i 




-1 




u{x,t) 



we obtain the corresponding operator A as 





"l 


0' 






-Ux 




6uUx —'iUxx 


A = -A 






d'x-Q 




dx - 







1 


Ux 






3uxx 6uux 



and the compatibility condition (4.3) yields the mKdV equation (2.5). 

V. The AKNS Method 

In 1973 Ablowitz, Kaup, Newell, and Scgur introduced [2,3] another method to determine an intcgrablc 
NPDE corresponding to a LODE. This method is now known as the AKNS method, and the basic idea 
behind it is the following. Given a linear operator X associated with the first-order system Vx = Xv, we are 
interested in finding an operator T (the operators T and X are said to form an AKNS pair) such that: 

(i) The spectral parameter A does not change in time, i.e. At = 0. 

(ii) The quantity Vt — Tv is also a solution to Vx = Xv, i.e. we have {vt — Tv)x = X{vt — Tv). 

(iii) The quantity Xt — Tx + — is a (matrix) multiplication operator, i.e. it is not a differential 
operator. 

Prom condition (ii) we get 

vtx - %v - Tvx = Xvt - XTv = {Xv)t - Xtv - XTv = {vx)t - Xtv - XTv = Vxt - XfV - XTv. (5.1) 

Using Vtx = Vxt and replacing Tvx by TXv on the left side and equating the left and right hand sides in 
(5.1), we obtain 

{Xt -Tx+XT- TX)v = 0, 

which in turn, because of (iii), implies 

Xt-Tx + XT -TX = 0. (5.2) 

We can view (5.2) as an integrable NPDE solvable with the help of the solutions to the direct and inverse 
scattering problems for the linear system Vx = Xv. Like (4.3), the compatibility condition (5.2) yields a 
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nonlinear evolution equation containing a first-order time derivative. Note that X contains the spectral 
parameter A, and hence T also depends on A as well. This is in contrast with the Lax method in the sense 

that the operator A does not contain A. 

Let us illustrate the AKNS method by deriving the KdV equation (2.1) from the Schrodinger equation 
(2.2). For this purpose we write the Schrodinger equation, by replacing the spectral parameter fc^ with A, 
as a first-order linear system Vx = Xv, where we have defined 



V :- 





, X 




"0 


u{x, t) 














1 









a 


/3" 






T = 










_P 


a 







"0 


0" 











Let us look for T in the form 



where the entries a, /3, p, and a may depend on x, t, and A. The compatibility condition (5.2) yields 

-ax - /3 + p{u - A) ut- Px + (^{u - A) - a{u - A) 
-Px + a-a ~(Jx + 13 - p{u- X) 

The (1, 1), (2, 1), and (2, 2)-entries in the matrix equation (5.3) imply 

P = -ax + {u- X)p, a = a-px, cTx = -ax- 
Then from the (1, 2)-entry in (5.3) we obtain 

ut + ^Pxxx - UxP - 2px{u - A) = 0. 



(5.3) 



(5.4) 



(5.5) 



Assuming a linear dependence of p on the spectral parameter and hence letting p — X( + p, in (5.5), we get 

2CxA^ + (^Cxxx - KxU + 2p,x - UxC \ X+iut + ^p.xxx - ^fJ-xU - UxM j = 0. 



Equating the coefficients of each power of A to zero, we have 



1 3 1 „ 

C = Ci, /i = -CiU + C2, Ut - -CiUUx - C2Ux + -ClUxxx = 0, 



(5.6) 



with ci and C2 denoting arbitrary constants. Choosing ci = 4 and C2 = 0, from (5.6) we obtain the KdV 
equation given in (2.1). Moreover, with the help of (5.4) we get 

a = Ux + C3, P = -4:X^ + 2Xu + 2u^ - Uxx, p = 4A + 2u, a = C3 - Ux, 
where C3 is an arbitrary constant. Choosing C3 = 0, we find 

Ux -4A^ + 2Au + 2u^ - Ux 
4A + 2u —Ux 



As for the Zakharov-Shabat system (2.4), writing it as Vx = Xv, where we have defined 

—iX u{x,t) 



X := 



—u{x,t) iX 
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we obtain the matrix operator T as 



T = 



and the compatibiUty condition (5.2) yields the NLS equation (2.3). 

As for the first-order linear system (2.6), by writing it as = Xv, where 

—iX u{x,t) 



X :-- 



-u{x,t) iX 



we obtain the matrix operator T as 



T 



-AiX^ + 2iXu'^ 



4A^u + 2iXux — Uxx 
UX^ - 2iAu2 



2u-^ 



-4A^u + 2iXux + Uxx + 2u^ 

and the compatibility condition (5.2) yields the mKdV equation (2.5) 
As for the first-order system Vx = Xv, where 

1 

-iX 

X := 



■Ux{x,t) 



1 / X 
^_-Ux{x,t) 



iX 



we obtain the matrix operator T as 



T = 



4A 



cos u sm u 
sin u — cos u 



Then, the compatibility condition (5.2) gives us the sine-Gordon equation 



Uxt = sm u. 



VI. Direct Scattering Problem 

The direct scattering problem consists of determining the scattering data when the potential is known. 
This problem is usually solved by obtaining certain specific solutions, known as the Jost solutions, to the 
relevant LODE. The appropriate scattering data can be constructed with the help of spatial asymptotics of 
the Jost solutions at infinity or from certain Wronskian relations among the Jost solutions. In this section 
we review the scattering data corresponding to the Schrodinger equation (2.2) and to the Zakharov-Shabat 
system (2.4). The scattering data sets for other LODEs can similarly be obtained. 

Consider (2.2) at fixed t by assuming that the potential u{x,t) belongs to the Faddcev class, i.e. u{x,t) 
is real valued and J^^dx{l + \x\) \u{x,t)\ is finite. The Schrodinger equation has two types of solutions; 
namely, scattering solutions and bound-state solutions. The scattering solutions are those that consist of 
linear combinations of e'*^^ and e"*''^ as a; ^ ±oo, and they occur for fc e R \ {0}, i.e. for real nonzero 
values of k. Two linearly independent scattering solutions fi and /r, known as the Jost solution from the left 
and from the right, respectively, are those solutions to (2.2) satisfying the respective asymptotic conditions 

fi{k,x,t) = e''''^ + o{l), f{{k,x,t) =ike'^'' + o{l), x ^ +oo, (6.1) 
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Mk, X, t) = e-"=^ + 0(1), f'Ak, X, t) = -ike-'''^ + o(l), 



-00, 



where the notation o(l) indicates the quantities that vanish. Writing their remaining spatial asymptotics in 
the form 



fi{k,x,t) 



fr{k, X, t) 



^ikx _^ L(fc,t)e-^*=^ 



T{k,t) 



+ 



T{k,t) 
R{k,t) 6**==^ 



+ 0(1), 



-oo. 



+00, 



(6.2) 



T{k,t) T{k,t) 

we obtain the scattering coefScients; namely, the transmission coefficient T and the reflection coefficients L 
and R, from the left and right, respectively. 

Let C+ denote the upper half complex plane. A bound-state solution to (2.2) is a solution that belongs 
to i-^(R) in the x variable. Note that i^(R) denotes the set of complex- valued functions whose absolute 
squares are integrable on the real line R. When u{x,t) is in the Faddeev class, it is known [5,7-9,16-19] 
that the number of bound states is finite, the multiplicity of each bound state is one, and the bound-state 
solutions can occur only at certain fc-values on the imaginary axis in C"*" . Let us use N to denote the number 
of bound states, and suppose that the bound states occur at fc = inj with the ordering < ki < • • • < kat. 
Each bound state corresponds to a pole of T in C"*". Any bound-state solution at k = iKj is a constant 
multiple of f\{iKj,x,t). The left and right bound-state norming constants c\j{t) and Crj{t), respectively, can 
be defined as 



cij{t) :-- 



f 

.J — ( 



dx f\{iK,j,x,tY 



^1/2 



f 

.J — ( 



-1/2 



and they are related to each other through the residues of T via 



Res (T, iKj) = i c\j {t) 7j {t) = i 
where the 7^ {t) are the dependency constants defined as 

fl{iKj,X,t) 
fii^lKj, X, t) 



Hit) 



(6.3) 



(6.4) 



The sign of 7j(t) is the same as that of (—1)''^ and hence Crj{t) = (—1)^ ^Ijit) cij(i)- 

The scattering matrix associated with (2.2) consists of the transmission coefficient T and the two 
reflection coefficients R and L, and it can be constructed from {KjljLi and one of the reflection coefficients. 
For example, if we start with the right reflection coefficient R{k, t) for A; e R, we get 



N 



r(fc,t)= m|±^| exp 



k 



{- f 



ds 



\og{l -\R{s,t)\^) 
s-k-iO+ 



fc e C+ u R, 



where the quantity iO+ indicates that the value for /c G R must be obtained as a limit from C"*". Then, the 
left reflection coefficient L[k,t) can be constructed via 



L{k,t) = - 



R{k,t)T{k,t) 



kGK. 



T(k,t) 

We will see in the next section that T{k,t) = T{k,0), \R{k,t)\ = |-R(fc,0)|, and \L{k,t)\ = \L{k,0)\. 
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For a detailed study of the direct scattering problem for the 1-D Schrodinger equation, we refer the reader 
to [5,7-9,16-19]. It is important to remember that u{x,t) for a; G R at each fixed t is uniquely determined 
[5,7-9,16-18] by the scattering data {i?, {k^}, {cy (t)}} or one of its equivalents. Letting Cj{t) := cy(t)-^, we 
will work with one such data set, namely {R, {kj}, {cj{t)}}, in Sections VII and VIII. 

Having described the scattering data associated with the Schrodinger equation, let us briefly describe 
the scattering data associated with the Zakharov-Shabat system (2.4). Assuming that u{x,t) for each t 
is integrable in x on R, the two Jost solutions a;,t) and (j){X,x,t), from the left and from the right, 
respectively, are those unique solutions to (2.4) satisfying the respective asymptotic conditions 



ip{X,x,t) 







0(1), 



+oo; 



(l){X,x,t) 







o{l), 



X — > — oo. 



(6.5) 



The transmission coefiicient T, the left reflection coefiicient L, and the right reflection coefficient R are 
obtained via the asymptotics 



V'(A,a;,i) 



L{\,t)e 



-iXx 



T{X,t) 



+ 0(1), 



-oo; (j){\x^t) = 



-iXx 



TiX,t) 
R{X,t) e*^^ 



o{l), 



+0O. (6.6) 



The bound-state solutions to (2.4) occur at those A values corresponding to the poles of T in C"*". Let us use 
{Xj}jLi to denote the set of such poles. It should be noted that such poles are not necessarily located on 
the positive imaginary axis. Furthermore, unlike the Schrodinger equation, the multiplicities of such poles 
may be greater than one. Let us assume that the pole has multiplicity Uj. Corresponding to the pole A^, 

one associates [4,20] nj bound-state norming constants Cjs(t) for s = 0,1, ... ,nj — 1. We assume that, for 
each fixed t, the potential u{x,t) in the Zakharov-Shabat system is uniquely determined by the scattering 
data {R, {Xj}, {cjs{t)}} and vice versa. 



VII. Time Evolution of the Scattering Data 

As the initial profile u{x, 0) evolves to u{x, t) while satisfying the NPDE, the corresponding initial 
scattering data 5(A, 0) evolves to S{X,t). Since the scattering data can be obtained from the Jost solutions 

to the associated LODE, in order to determine the time evolution of the scattering data, we can analyze the 
time evolution of the Jost solutions with the help of the Lax method or the AKNS method. 

Let us illustrate how to determine the time evolution of the scattering data in the Schrodinger equation 
with the help of the Lax method. As indicated in Section IV, the spectral parameter k and hence also the 
values Kj related to the bound states remain unchanged in time. Let us obtain the time evolution of f\{k, x, t), 
the Jost solution from the left. Prom condition (ii) in Section IV, we see that the quantity dtf\ — Af\ remains 
a solution to (2.2) and hence we can write it as a linear combination of the two linearly independent Jost 
solutions /i and /r as 

dth{k, X, t) - {-Adl + Qud^ + 3u^) /i(fc, x, t) = p{k, t) fi{k, X, t) + q{k, t) /,(fc, x, t), (7.1) 

where the coefficients p{k,t) and q{k,t) are yet to be determined and A is the operator in (4.7). For each 
fixed t, assuming u{x, t) = o(l) and Ux{x, t) = o(l) as a; — > +oo and using (6.1) and (6.2) in (7.1) as a; — > +oo, 
we get 

ate*'^^ + 4a3e*'=^ = p{k, t) e*'=^ + q{k, t) 



T(k,t) 



ikx _|_ -^(^' ^) gifex 



T(k,t) 



+ 0(1). 



(7.2) 
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Comparing the coefficients of e'*^^ and e ''^^ on the two sides of (7.2), we obtain 

q{k,t) = 0, p{k,t) = -Aik^. 
Thus, fi{k,x,t) evolves in time by obeying the hnear third-order PDE 

dth-Ah = -Ueh. (7.3) 
Proceeding in a similar manner, we find that /r(A;, x, t) evolves in time according to 

dtf. - Af, = Aik^fr. (7.4) 

Notice that the time evolution of each Jost solution is fairly complicated. We will see, however, that the 
time evolution of the scattering data is very simple. Letting x — oo in (7.3), using (6.2) and u{x, t) = o(l) 
and Ux{x,t) = o(l) as a; ^ — oo, and comparing the coefficients of e^''^ and e~^''^ on both sides, we obtain 

dtT{k, t) = 0, dtL{k, t) = -8ik^L{k, t), 

yielding 

T(fc, t) = T(k, 0), L{k, t) = L{k, 0) e"^**^'*. 
In a similar way, from (7.4) as a; ^ +oo, we get 

R{k,t) = R{k,Q)e^'^^K (7.5) 

Thus, the transmission coefficient remains unchanged and only the phases of the reflection coefflcients change 
as time progresses. 

Let us also evaluate the time evolution of the dependency constants ^j{t) defined in (6.4). Evaluating 
(7.3) at fc = iKj and replacing f\{iKj,x,t) by 'yj{t)fi{iKj,x,t), we get 

fr{inj,x,t)daj{t) + 'jj{t)dtfr{iKj,x,t) - jj{t)Afr{iKj,x,t) = -4:K^j'yj{t)f,{iK,j,x,t). (7.6) 

On the other hand, evaluating (7.4) at fc = inj , we obtain 

■lj{'t)9tfr{iK.j,X,t) - Jj{t)AfriiKj,X,t) = 4:K^jJj{t) fr{iK,j , X,t). (7.7) 

Comparing (7.6) and (7.7) we see that dtjj{t) = — 8K^7j(t), or equivalently 

7,W=7,(0)e-««?*. (7.8) 
Then, with the help of (6.3) and (7.8), we determine the time evolutions of the norming constants as 

cy it) = cy (0) e^'^?*, Cr, it) = Cr, (0) e'^^'*. 

The norming constants Cj{t) appearing in the Marchenko kernel (8.1) are related to c\j{t) as Cj{t) := cij{t)^, 
and hence their time evolution is described as 

c,(t) = c,(0)e«'^?*. (7.9) 
11 



As for the NLS equation and other integrable NPDEs, the time evohition of the; related sc;attc;ring data 
sets can be obtained in a similar way. For the former, in terms of the operator A in (4.8), the Jest solutions 
ip{\,x,t) and (j){X,x,t) appearing in (6.5) evolve according to the respective linear PDEs 

ijJ^-A^ = -2i\^^, (j)t-A(j) = 2iX'^(j). 

The scattering coefBcients appearing in (6.6) evolve according to 

T(A,i) =T(A,0), i?(A,t) =i?(A,0)e^'^'*, L{\,t) = L{X,0)e-'^'^^*. 



(7.10) 



Associated with the bound-state pole Xj of T, we have the bound-state norming constants cjs {t) appearing 
in the Marchenko kernel fl{y,t) given in (8.4). Their time evolution is governed [4] by 



[Cj(«,-l)(i) 9("j-2)W ••• Cjo{t)] = [Cj(n,-l)iO) Cj(„^._2)(0) 

where the rij x rij matrix Aj appearing in the exponent is defined as 



Cjo(0)]e 



(7.11) 



Ar.= 



-iXj 









-iXj -1 



-1 
-iX 




-iXi 









—iXi 



1 





-1 



-iXi 



VIII. Inverse Scattering Problem 

In Section VI we have seen how the initial scattering data ^(A, 0) can be constructed from the initial 
profile u{x,0) of the potential by solving the direct scattering problem for the relevant LODE. Then, in 
Section VII we have seen how to obtain the time-evolved scattering data S{X,t) from the initial scattering 
data S{X,0). As the final step in the 1ST, in this section we outline how to obtain u{x,t) from S{X,t) by 
solving the relevant inverse scattering problem. Such an inverse scattering problem may be solved by the 
Marchenko method [5,7-9,16-19]. Unfortunately, in the literature many researchers refer to this method as 
the Gel'fand-Levitan method or the Gel'fand-Levitan-Marchenko method, both of which are misnomers. The 
Gel'fand-Levitan method [5,7,16,17,19] is a different method to solve the inverse scattering problem, and the 
corresponding Gel'fand-Levitan integral equation involves an integration on the finite interval (0, x) and its 
kernel is related to the Fourier transform of the spectral measure associated with the LODE. On the other 
hand, the Marchenko integral equation involves an integration on the semi-infinite interval (x,-|-oo), and its 
kernel is related to the Fourier transform of the scattering data. 

In this section we first outline the recovery of the solution u{x,t) to the KdV equation from the cor- 
responding time-evolved scattering data {R,{Kj},{cj{t)}} appearing in (7.5) and (7.9). Later, we will 
also outline the recovery of the solution u{x, t) to the NLS equation from the corresponding time-evolved 
scattering data {R, {Xj} , {cjs{t)}} appearing in (7.10) and (7.11). 

The solution u{x,t) to the KdV equation (2.1) can be obtained from the time-evolved scattering data 

by using the Marchenko method as follows: 

(a) From the scattering data {R{k, t),{Kj}, {cj{t)}} appearing in (7.5) and (7.9), form the Marchenko kernel 
O defined via 

^iy,t)-=7r dkR{k,t)e'''y + y^c,{t)e-^^y. (8.1) 
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(b) Solve the corresponding Marchenko integral equation 

K{x,y,t) + n{x + y,t)+ dz K{x,z,t)Q.{z + y,t) = Q, x<y<+oo, (8.2) 



and obtain its solution K{x, y, t). 
(c) Recover u{x, t) by using 



u{x,t) 



dK{x, X, t) 
dx 



(8.3) 



The solution u{x, t) to the NLS equation (2.3) can be obtained from the time-evolved scattering data 
by using the Marchenko method as follows: 

(i) Prom the scattering data {i?(A, t), {Aj}, {cjs(f)}} appearing in (7.10) and (7.11), form the Marchenko 
kernel O as 

Q{y,t) := - / dXR{X,t)e'>^y + Y^ ^ c,-.(t)^ e'^^''. (8.4) 

(ii) Solve the Marchenko integral equation 

K{x,y,t) — ^{x + y,t) + / dz ds K{x,s,t)il{s + z,t)Q{z + y,t) = 0, x<y<+oo, 

J X J X 

and obtain its solution K{x,y,t). 

(iii) Recover u{x,t) from the solution K{x,y,t) to the Marchenko equation via 

u{x,t) = —2K{x,x,t). 

(iv) Having determined K{x,y,t), one can alternatively get |M(a;,f)p from 

dG{x, x, t) 



\u{x,t)y 



dx 



where we have defined 



G{x 



/•oo 

,y,t):=— dz K{x, z,t)^l{z + 

J X 



y,t). 



IX. Solitons 

A soliton solution to an integrable NPDE is a solution u{x, t) for which the reflection coefficient in the 

corresponding scattering data is zero. In other words, a soliton solution u{x,t) to an integrable NPDE is 
nothing but a rcflcctionlcss potential in the associated LODE. When the reflection coefficient is zero, the 
kernel of the relevant Marchenko integral equation becomes separable. An integral equation with a separable 
kernel can be solved explicitly by transforming that linear equation into a system of linear algebraic equations. 
In that case, we get exact solutions to the integrable NPDE, which are known as soliton solutions. 

For the KdV equation the A''-soliton solution is obtained by using R{k,t) = in (8.1). In that case, 
letting 

ci{t) e-^^y 



X{x) := [e-'^i^ e- 



Y{y,t) :-- 



C2{t)i 
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we get fl{x + y,t) = X{x)Y{y,t). As a result of this separability the Marchenko integral equation can be 
solved algebraically and the solution has the form K{x, y, t) = H{x, t) Y{y, t), where H{x, t) is a row vector 
with A'' entries that are functions of x and t. A substitution in (8.2) yields 



K{x,y,t) = -X{x)r{x,t)-'Y{y,t), 
where the N x N matrix T{x,t) is given by 

T{x,t):=I+ dzY{z,t)X{z), 

J X 

with I denoting the N x N identity matrix. Equivalently, the {j, /)-entry of T is given by 



(9.1) 



(9.2) 



Tji = 6ji + 



Kj + Kl 



with 6ji denoting the Kronecker delta. Using (9.1) in (8.3) we obtain 



u{x, t)=2— [X{x) T{x, t)-^Y{x, t)] = 2 tr — [Y{x, t) X{x) T{x, t)-^] , 



dx 



where tr denotes the matrix trace (the sum of diagonal entries in a square matrix). Prom (9.2) we see that 
—Y{x, t) X{x) is equal to the a;-derivative of r(x, t) and hence the A'-soliton solution can also be written as 



u{x,t) 



-2tr — 
ox 



dT{x,t) 
dx 



nx,t)-' 



dx 



^detT{x,t) 



dx 



det r(a;, t) 



(9.3) 



where det denotes the matrix determinant. When A' = 1, we can express the one-soliton solution u{x,t) to 
the KdV equation in the equivalent form 

u{x, t) = -2 kI sech^ {kiX - 4K;f f + 6*) , 



with e := log i/2ki/ci(0). 

Let us mention that, using matrix exponentials, we can express [6] the A'-soliton solution appearing in 
(9.3) in various other equivalent forms such as 



where 



u{x,t) = -ACe-^'=+^^'*T{x,t)-^AT{x,t)-^e-'^'=B, 

A := diag{Ki, K2, • • • , kn}, 
B^:=[l 1 ... 1], C:=[ci{0) C2(0) ... cjv(O)]. 



(9.4) 



Note that a dagger is used for the matrix adjoint (transpose and complex conjugate), and B has N entries. 
In this notation we can express (9.2) as 



T{x,t) 



f 

J X 



dze-^^BCe-'^e^'^\ 



As for the NLS equation, the well-known A'-soliton solution (with simple bound-state poles) is obtained 
by choosing i?(A. f.) = Q and Uj = 1 in (8.4). Proceeding as in the KdV case, we obtain the A'-soliton solution 
in terms of the triplet A, B, C with 



A := diag{-iAi, -iX2, . . . , -iXn}, 



(9.5) 
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where the complex constants Xj are the distinct poles of the transmission coefBcient in C"*", B and C are 
as in (9.4) except for the fact that the constants Cj(0) are now allowed to be nonzero complex numbers. In 
terms of the matrices P{x,t), M, and Q defined as 



P{x, t) := diag{e2^^i^+4'^?*, e^'^^x+^A^t^ _ _ _ ^ ^2i\MX+Ai\lt^^ j^,^ .__ 



Xj - Xi 



> Qji := 



-ICjCl 

Xj - Xi 



we construct the iV-soliton solution u{x, t) to the NLS equation as 

u{x,t) = -2Bt [I + P{x,tyQP{x,t) M]~'^ P{x,t)^C^ , 

or equivalently as 



i{x,t) = -2Bte-^t-r(a;,i)-ie-^'^+4^(^')'*Ct, 



where we have defined 
T{x,t) :=! + 



Using (9.4) and (9.5) in (9.8), we get the (j, /)-entry of r(x,t) as 



(9.6) 
(9.7) 

(9.8) 



m=l i^m - Xj){Xm - Xi) 



Note that the absolute square of u{x, t) is given by 

\u{x,t)\^ = tr 



d_ 
dx 



_d_ 

dx 



detr(a;,t) 



det r(a;, t) 

For the NLS equation, when N = 1, from (9.6) or (9.7) we obtain the single-soliton solution 

-8ci(Im[Ai])2 e-2iAlx-4»(Ai)^t 



u{x,t) = 



4(Im[Ai])2 + |ci|2 e-4x(Im[Ai])-8t(Im[A?]) 



where Im denotes the imaginary part. 
X. Future Directions 

There are many issues related to the 1ST and solitons that cannot be discussed in such a short review. 
We will briefly mention only a few. 

Can we characterize integrable NPDEs? In other words, can we find a set of necessary and sufficient 
conditions that guarantee that an IVP for a NPDE is solvable via an 1ST? Integrable NPDEs seem to 
have some common characteristic features [1] such as possessing Lax pairs, AKNS pairs, soliton solutions, 

infinite number of conserved quantities, a Hamiltonian formalism, the Painleve property, and the Backlund 
transformation. Yet, there docs not seem to be a satisfactory solution to their characterization problem. 

Another interesting question is the determination of the LODE associated with an 1ST. In other words, 
given an integrable NPDE, can we determine the corresponding LODE? There does not yet seem to be a 
completely satisfactory answer to this question. 

When the initial scattering coefficients are rational functions of the spectral parameter, representing the 
time-evolved scattering data in terms of matrix exponentials results in the separability of the kernel of the 
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Marchenko integral equation. In that case, one obtains explicit formulas [4,6] for exact solutions to some 
integrable NPDEs and such solutions are constructed in terms of a triplet of constant matrices A, B, C 
whose sizes are p x p, p x 1, and 1 x p, respectively, for any positive integer p. Some special cases of such 
solutions have been mentioned in Section IX, and it would be interesting to determine if such exact solutions 
can be constructed also when p becomes infinite. 
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